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ABSTRACT 

By introducing the concept of the age distribution of cosmic 
rays it is possible to decouple the spatial from the momentum 
transport; and simple leaky-box type equations result. The 
influence of spatial inhomogeneities, geometries and source 
distributions enters the spatially homogeneous, infinite 
(i.e. leaky box) problem through appropriate mean lifetimes. 

A precise prescription of how to obtain these mean lifetimes: 
for comparison with data measured in the vicinity of the 
solar system they have to be calculated from the age 
distribution at the solar system. 

1. Introduction. The propagation of cosmic ray particles in the Galaxy 
at large momenta, p > 10 GeV/c, where spatial diffusion in partially 
random magnetic fields dominates convection and adiabatic deceleration in 
the galactic wind, is described by the steady-state transport equation 
(e.g. Schlickeiser 1983) 


V • [K(r,p) Vf] + L f + q(r) Q(p) = 0 ( 1 ) s 

where f is the particle ' s phase space density, K(r,p) = K(p) K(r) the 
spatial diffusion tensor, with K(p) chosen dimensionless without loss of 
generality, and Lp is the momentum operator 


-2 d 


[ P 2 D(p) •— - 


describing first (pgain) an d second (D ( p ) ) order Fermi acceleration as 
well as continuous Tp^ oss ) and catastrophic (x„) momentum loss processes. 
q(j) Q(p) represents sources and sinks of particles. 


2. Decoupling of spatial and momentum transport. If the parameters 
D(p), f'gaj.n, floss an< d T c do n °t- depend on position (r) , equation (1) is 
solved by the convolution (Lerche and Schlickeiser 1985a) 

CO 

f(p,r) = / dy G(p,t) N(r,y) ( 3 ) } 

0 

with 



G(p,y) 


3G(p,y) 

3y 


0 


( 4 ) 


and G(p,y) satisfying given momentum boundary conditions and 
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(5) 


Lim G(p,y) = 0 

y-KX) 


G(p,y=0) = Q(p)/K(p) 

(6) . 

N(r,y) satisfies given spatial boundary conditions, and 


7 • [K(r) V N(r,y)J = 

~ ~ 9 y 

(7), 

N(r ,y = 0) = q(r) 

(8), 

and N(r,y) is normalized as 


oo 

/ dy N(r,y) = 1 

(9). 


0 

N(r,y) is coramonly referred to as the age distribution of cosmic rays. 


The solution of the system of equations (4) “ (8) would yield via (3) 
the particle’s momentum spectrum f(p,r) at all positions r in the Galaxy. 
However, particle flux measurements can only be obtained at the position 
of the solar system. In order to compare with these observations we 
would infer from the general solution f(p,r) the momentum spectrum at the 
solar system, 

oo 

H(p) = f(p,r ) = | dy G(p,y) N(r g ,y) (10) 

0 

by setting r = r g in (3). 


3. Leaky-box equations. However, instead of following this scheme to 
solve the two partial differential equations (4) and (7), it is much 
easier to derive from (4) a simpler equation for the momentum spectrum 
H(p) at the solar system from the known cosmic ray age distribution 
N(r g ,y) at the solar system. Integrate equation (4) after multiplying it 
with N(r s ,y) to obtain 


L p / dy G(p,y) N(r g ,y) = jf dy - G 4j ^ - N(r g ,y) 


Q(p) N(r ,y=0) « 9N(r ,y) 

~ S - - J dy G(p,y) 


K(p) 


0 


9y 


( 11 ), 


where the right hand side has been integrated by parts taking into ac- 
count (5) and (6). With (10) equation (11) reads 


OO 

L H(p)+K(p) J dy G(p,y) 
p O 


9N(r s »y^ 

9y 


- Q(p) N(r ,y=0) 


which is an ordinary differential equation in momentum. 


( 12 ) 


If the age distribution at the solar system is a pure exponential, 
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N( Hs’ y) = <r >(r~) ~ ex P (-y/<T>(r a » (13), 

one infers N(r s ,y=0) = 1/<T>(r s ) and . 3N(r s ,y) /9y = (-1/<T>(r s )) N(r s ,y), 
so that (12) attains the form of a simple leaky— box equation 


L p H(p) 


K(p) H(p) _ Q(p) 

<T>(r s ) “ <T> (r s ) 


(14), 


where the whole spatial problem enters through the constant parameter 
< T > (r s ) characterizing the exponential age distribution (13). It has to 
be emphasized that equation (14) is only valid at the position of the 
solar system. If for some reason, one wants to calculate the momentum 
spectrum of particles at some other place in the Galaxy, say r^, one has 
to calculate new parameters <T>(r£) from the respective age distribution 
N(r i>y ). 

But is the age distribution of cosmic rays arriving at the solar 
system pure exponential, as assumed in (13) and (14)? As shown by Lerche 
and Schlickeiser (1985b), in a variety of circumstances the age distri- 
bution at the solar system can be expanded in a series of exponentials: 

N( 5 S ’ y) = l <t 7~ ex P (-y/< T > m ) (is). 

m=1 m m 


The convolution of N(r s ,y) from (15) with equation (4) yields for 
particle momentum spectrum near the solar system 

00 

H(p) = y c H (p) 

L , m m r 
m=1 


(16), 


where H m (p) is the solution of equation (14) with <T> m playing the role 
of <T>(r s ). So in general, H(p) is attained by an (infinite) sum of 
leaky-box solutions. Very often it suffices to consider only the m = 1 
contribution, when the first eigenvalue dominates the behaviour of the 
age distribution. 


4. Conclusions. In studies of the transport and acceleration of cosmic 
energetic particles the age distribution plays the same role as the 
escape probability in studies of photon propagation through hot plasma 
(see e.g. Sunyaev and Titarchuk 1980). The influence of spatial in- 
homogeneities, geometries and source distributions can be estimated using 
the solution for the spatially homogeneous infinite problem (leaky-box 
equation (14)) with the appropriate escape lifetimes <T>. A precise 
prescription of how to obtain these escape lifetimes has been given: 
namely, they have to be calculated from the age distribution at the 
spatial position of the observer. 
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